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Introduction 

The radial velocity method for detecting extra-solar planets relies on measuring the star's wobble around the 

system's center of mass. Since this is an indirect method, we may ask if there are other dynamical effects that 

can mimic such wobble. In recent articles|T| |2] |3], we modeled the effect of a nearby binary system on a star's 

("^ ■ radial velocity. We showed that, if we are unaware of this nearby binary, for instance because one component is 

CS| I unresolved or both components are faint stars, the binary's effect may mimic a planet. Here, we review this work, 

^ i explaining in which circumstances the binary's effect may mimic a planet and we discuss what can be done in 

^ ' practice in order to distinguish between these two scenarios (planet or nearby binary). 

CN I 1. Overview of theory and results 

n , I We study a triple system composed by a star with mass m2, at distance r2 from the centre of mass of a binary 

nQ ' with masses mo + mi, and inter-binary distance ri. We assume that p = |ri|/|r2| ^ 1 (hierarchical system). The 

_; I Hamiltonian in Jacobi coordinates is |I3| 

Ph. Gmomi G(mo + mi)m2 

I , rt = \- r , (I) 
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where an approximation of the perturbation term is 



Gm2 momi p^ ^ 

F = ■ -(Scos-'^-l) , (2) 

fl mo + mi 2 

^ . with S the angle between rj and r2- 

OO ' From Eq. ([U we see that the motion is, approximately, a composition of two Keplerian orbits described by: ri, 

with semi-major axis ai, eccentricity ei, and period Ti =2 n/n\\ and r2, with semi-major axis a\, eccentricity ey, 
^^ ' and period T\=2 n/ni. 

The radial velocity of the star m2 is Vr — Vrk + Vrp, where Vrk is a Keplerian term that describes the motion 

around a "star", of mass m^ + mi, located at the binary system's centre of mass, and Vrp is a small perturbation||T]|2l- 

^^ . The radial velocity data of the star m2 is first fitted with a Keplerian radial velocity curve, Vrk- After subtracting 

Vrk, we are left with the perturbation term, Vrp, from which we could in principle infer the presence of the nearby 

binary. However, as we will show next, this is not always possible in practice. 



1 . 1 Coplanar circular orbits 

In the case of coplanar circular orbits we havelH] 

Vrp = Kocos{n2t + eo)+Kicos{{2ni-3n2)t + ei)+K2Cos{{2ni-n2}t + e2) ■ (3) 

The term with frequency n2 is incorporated in the main Keplerian curve, Vrk- Since ni ^n2, the 2nd and 3rd terms 
have very close frequencies that can only be resolved if the observation timespan t,yhs > ^2/2. If there is enough 
resolution and precision, we identify both signals and we conclude that they should not be caused by planets (as 
such close orbits would be unstable). However, since j/Ti | = 5 \K2\, in practice, due to limited precision, we may 
only be able to observe the signal with frequency 2 nj — 3 n2- In this case, we may think there is a planet companion 
to the observed star. 

We simulated a coplanar triple system composed of m2 = Mq on a circular orbit of period T2 = 22 y, around 
a binary mo = 0.7 Mq and mi ~ O.35M0 with a circular orbit of period Ti =411 d. We computed radial velocity 
data points over ?„/„ = 11 y at precision 0.8 m/s. In Fig. 1 we see a signal with frequency 2ni — 3«2 (amplitude 
0.9 m/s at 223 d) that mimics a planet of ISAfe. 




Figure 1: Periodogram of residuals leftover after removing 
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Figure 2: Periodogram of residuals leftover after removing Vrk- 



1 .2 Non-coplanar circular orbits 

In the case of non-coplanar circular orbits, Vrp is a combination of 6 periodic terms|l2] with frequencies: ni, 
3ni, 2ni ±n2 and 2ni ±3n2- Depending on the observation's precision and resolution we may observe one or 
more of these terms. If all these have well separated frequencies we mistake them by planet(s). 

We simulated a triple system composed of m2 — Mp, on a circular orbit of period T2 — 4.2 y, around a binary 
niQ — mi = O.25M0 with a circular orbit of period T\ — 85 d. The relative inclination is / — 30°. We computed 
radial velocity data points over tohs = 1 1 y at precision 0.7 m/s. In Fig. 2 we see signals with frequency 2«i — 3 n2 
(amplitude 0.8 m/s at 46 d), and frequency 3 ^2 (amplitude 1 .4 m/s at 5 16 d) that mimic planets oil Me and 20M£. 



1 . 3 Eccentric coplanar orbits 



It can be shown that, generally, Vrp is a composition of short period terms||2]. When the orbits are circular we 
saw above that there are finite number of periodic terms. When the orbits are eccentric we must express Vrp as 
an expansion in ei and 62, hence there is an infinite number of periodic terms. However, in practice only a finite 
number of periodic terms (the ones that appear at lowest order in ei and 62) are important. In the eccentric 2D 
case, up to 1st order in the eccentricities, there are 12 frequencies[|2|: n2, 2ni —n2, 2ni -~3n2, ni ~n2, n\ — 3«2, 
3ni — «2> 3«i — 3«2i «i +"2^ 2n2, 2ni, 2n\ —An2 and 2n\ — 2n2- 

We simulated a coplanar triple system composed of a star OT2 = M© with 62 =0.1 and period T2 = 22 y, in 
orbit around a binary with mq = 0.7 Mq and mi = O.35M0, ei = 0.2 and period Ji =411 d. We computed radial 
velocity data points over ?„/„ = 11 y at precision 0.8 m/s. In Fig. 3 we see signals with frequencies 2ni — 3n2 
(amplitude 0.8 m/s at 223 d) and rti — 3n2 (amplitude 1.4 m/s at 487 d). Since n2 <C mi these mimic planets of 
15 Me and 3AMe at the 2/1 mean motion resonance. 




Figure 3: Periodogram of residuals leftover after removing Vrk- 



1 .4 Secular evolution 

If the star 1112 has an eccentric orbit (e2 7^ 0) and if ?„/„ > T2, signals at or nearby harmonics of «2 appear^ E] . 
These may be mistaken by planet(s) near a mean motion resonance with a companion "star" of mass mo + miQ. 
However, as t^bs increases, the short period terms should become negligible with respect to the orbits' secular 
evolution. 

The secular evolution can be obtained by averaging Eq. ([T]i with respect to the orbital periods Ti and 72. 
Moreover, since ai -C 02 (hierarchical system) and m2 > mo + mi (perturbing binary less massive than observed 
star), it can be shown that the star's motion around the binary's center of mass coincides, approximately, with the 
invariant plane [ 3 1 . The secular Hamiltonian, obtaining after averaging equation (2) over the mean motion of both 
orbits, is approximateljo 



G niQini 



m2 



16mf) + mi (1 —e^y/^ , 



[(2 + 3ei)(3cos2/-l) + 15eisin2/cos(2«i)] 



(4) 



with ei (binary's eccentricity), coi (binary's argument of pericentre), / (relative inclination). 

From Hamilton's equations we obtain the secular evolution of the star's orbit. Since Eq. (|4|l does not depend 
on CO2 (star's argument of pericentre) then 62 = and 



3 , mnini 

(02 = -A- -:r 

4 (mo + mi)'^ 



"2 



(l-.2)2' 



where 



. ,1 3 
A= - + -, 

' 2 4 
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•1 I (3cos^/-l)H ei(l-cos2/)cos(2a)i) 



(5) 



(6) 



depends on the binary's (mo + mi) orbit. 

In Fig. 4 we show A calculated at fixed eii3|. If / < 40° we can eliminate C0\ from Eq. (|6]l since (Oi circulates. 
If / > 40° we place the binary at the Kozai equilibrium (ei = -^/l — (5/3) cos^/, Oi = ±90°). 

Therefore, Vr is better fitted by the radial velocity curve of a precessing ellipse 13] 



Vr = 



n2a2 nih sin/2 
(1 — gj)^^^ m2+mh 



(cOs(/2 + 0)2,0 + Oht) + e2 COS(C02,0 + 0>2t)) 



(7) 



where m^ = mo + mi and 6)2 is given by Eq. (|5]). 

We simulated a system composed of a star m2 = Mq with an eccentric orbit (e2 = 0.2) of period T2 = 4.24 y, 
and a binary with masses mo = O.35M0 and mi = O.15M0, period Ti = 85 d, relative inclination / = 60° and 
eccentricity ei = 0.76 (at the Kozai equilibrium). We computed radial velocity data points over f„ij = 6 y at 
precision 5 m/s. In Fig. 5 (top) we see the residuals leftover after fitting a fixed Keplerian radial velocity curve. 
There is an an obvious peak at 606 day which is a frequency close to an harmonic of n2- In Fig. 5 (low) we see the 
residuals leftover after fitting a radial velocity curve with precession (Eqs. l7l5l l. We obtain (02 = — 0.25°/y and the 
peak at 606 d is no longer important. 



'in (3) we present a longer but more correct derivation of the triple system's secular evolution. 
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Figure 4: The normalized precession rate (A). 




Figure 5: Periodogram of residuals leftover after fitting Keplerian radial velocity curve (top) or radial velocity 

curve of precessing ellipse (low) to the data. 



Conclusion 

When f„ftj < r2, we observe short-period terms due to the binary that may be mistaken by planet(s). In order to 
distinguish the binary's effect from planet(s) we need precise observations over a reasonable long timespan, which 
is often not possible in practice. However, a signal with frequency, «p/, and amplitude, Kpi, mimics a planet with 
parameters Uj2J 

P' ~ 2/3 ' ^ ' 

"pl 

. J „ mo + mi +1112 ,o^ 

mpi sm I pi = Kpi . (9) 

ripi a pi 

Therefore, when the detected signals' frequencies are well separated we can invert these expressions to predict 
binary system's parameters and check if they are realistic [Tl|2l. 

When fofov :§> T2 the orbits' secular evolution dominates over the short-period terms. In this case, we showed that 
the star's orbit precesses due to the binary, and measuring this precession rate provides estimates for the binary's 
parametersl3|. 

We detected such a precessional effect in the radial velocity data of the star v-Octantis A[f3l. This could 
provide an alternative explanation for the signal at 417 d that had previously been identified as a planetQ. Since 
V-Octantis A has a nearby companion (v-Octantis B) on a 2.8 y orbitf4l, a planet at 417 d cannot be stable on 
a coplanar prograde orbit (although it can survive on a coplanar retrograde orbit)[5J. We measured retrograde 
precession of — 0.86°/y which could be explained if v-Octantis B was a double star inclined / > 45° with respect to 
the main binary's orbitfj]. Moreover, we observed that after fitting a precessing orbit to the radial velocity data, the 
signal at 417 d was no longer prominent in the leftover residuals [3j. Although, retrograde precession of the main 
binary could also be explained by a planet on a highly inclined orbit||3l, these orbits do not seem to be stable||3]. 
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